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Understanding the performance of the low energy neutrino factory: the dependence 

on baseline distance and stored-muon energy 

Peter Ballett 1 B and Silvia Pascoli^Q 
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We study the physics reach of a Low Energy Neutrino Factory (LENF) and its dependence on the 
chosen baseline distance, L, and stored-muon energy, E^, in order to ascertain the configuration of 
the optimal LENF. In particular, we study the performance of the LENF over a range of baseline 
distances from 1000 km to 4000 km and stored-muon energies from 4 GeV to 25 GeV, connecting 
the early studies of the LENF (1300 km, 4.5 GeV) to those of the conventional, high-energy neutrino 
factory design (4000 km and 7000 km, 25 GeV). Three different magnetized detector options are 
considered: a Totally- Active Scintillator Detector (TASD) and two models of a liquid-argon detec- 
tor distinguished by optimistic and conservative performance estimates. In order to compare the 
sensitivity of each set-up, we compute the full ^-dependent discovery contours for the determination 
of #i3 7^ 0, Sep ^ {0, it} and sign(Amf 3 ). For large values of #13, as recently confirmed by the Daya 
Bay and RENO experiments, the LENF provides a strong discovery potential over the majority of 
the L-E^ parameter space and is a promising candidate for the future generation of long baseline 
experiments aimed at discovering CP-violation and the mass hierarchy, and at making a precise 
determination of the oscillation parameters. 

PACS numbers: 14.60. Pq 



I. INTRODUCTION 

Over the last 50 years, neutrino oscillations have ad- 
vanced from a tentative hypothesis to a well confirmed 
physical phenomenon and have been observed in vari- 
ous channels and over a wide range of energy scales. 
Until recently, the results of several experiments (see, 
e.g. [l[) had only been able to provide us with measure- 
ments of four of the six oscillation parameters: 612, #23, 
Am^ and lAmjj-J [2|-(5j]. This situation has changed over 
the last few months, with new experimental searches 0- 
succeeding in measuring the last mixing angle #13 
and confirming that it is large with central values re- 
ported by Daya Bay and RENO of sin 2 26 13 = 0.092 and 
sin 2 20i3 = 0.113, respectively [IH H(i||. The discovery of 
non-zero #13 has increased the likelihood that the res- 
olution of the final uncertainties in the neutrino mixing 
parameters will be possible within a medium-term exper- 
imental program. There are two particularly important 
unknowns that will dominate this project: measuring the 
complex phase S, which governs CP-violating effects in 
the leptonic sector and is currently unconstrained, and 
determining the sign of Am 2 1; which dictates the type of 
neutrino mass-ordering and controls the enhancement of 
the oscillation probabilities in the Earth due to matter 
effects for long-baseline and atmospheric neutrinos. 

The values of the unknown neutrino mixing parame- 
ters may well turn out to have decisive implications for 
potential extensions of the Standard Model and for un- 
derstanding the problem of flavor. Improving our knowl- 



edge of these parameters is one of the main aims of 
the next generation of long-baseline oscillation searches. 
These experiments are primarily aimed at studying the 
sub-dominant oscillation probability P{v e — > v^), or 



These probabilities have a complex depen- 
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dence on the oscillation parameters and their extraction 
from data suffers from the problem of degeneracies [TiT - 
[l9j . Various strategies have been devised to weaken or 
resolve this problem [2CH22| . for example the energy de- 
pendence of the degenerate solutions can be exploited 
using a wide- band beam [23|, [24[ , matter effects can be 
used to lift the hierarchy degeneracy [III [26| whilst, for 
the measurement of 623 and its deviations from maximal- 
ly, the low-energy part of the spectrum is particularly 
important [2?5 as well as a precise determination of the 
disappearance probability. 

Amongst future long-baseline neutrino experiments, 
the Neutrino Factory (NF) [28|-|30| is a leading design 
with excellent physics reach. At a NF, muons (or anti- 
muons) produced via pion decay are accelerated to a com- 
mon energy (originally around 25 GeV) and are in- 
jected into a storage ring. A significant number of high- 
gamma muons will decay producing a highly collimated 
beam of neutrinos, whose spectral content is known with 
high precision. The decay of yT (p + ) in the storage 
ring produces an initial beam with two neutrino com- 
ponents: and v e (F M and v e ). At the detector, two 
different signals will be present: for instance with /i~ de- 
cays, the right-sign muon events which derive from the 
observation of coming from the disappearance chan- 
nel and the wrong-sign muon events which are due to 
v& oscillations. It is necessary to distinguish these 

two signals in order to reconstruct the appearance prob- 
ability. This requires a magnetized detector which is able 
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to distinguish fi + - and u~-events. The baseline choice for 
the conventional NF [3l[ is a Magnetized Iron Neutrino 
Detector (MIND) which provides excellent background 
rejection and very good energy resolution but low de- 
tection efficiency for neutrinos with energies in the few 
GeV range [H, |33| . Other options include a magnetized 
Totally- Active Scintillator Detector (TASD) [30j or a de- 
tector based on the liquid-argon time projection cham- 
ber [33 - f37| . The latter detectors provide excellent effi- 
ciency for neutrinos with energies as low as a few hun- 
dred MeV to 1 GeV, excellent energy resolution and low 
backgrounds but their magnetization is extremely chal- 
lenging nn. 

Before the measurement of 6*13, the baseline configura- 
tion of the NF [3l| used muons with an energy of 25 GeV 
and two different baselines, at approximately 4000 and 
7500 km, with two MIND detectors, a 100 kton one at 
the shorter baseline and a 50 kton one at the 'magic' 
baseline [H, Ho| . The combination of two distinct base- 
line distances was designed to help resolve the problem of 
parameter degeneracies by providing complementary in- 
formation of the oscillation probability at different points 
in parameter space. This set-up was optimized for small 
values of #13 and has been shown to have an excellent 
physics reach to #13, CP-violation and the mass hierar- 

chy PEHQ3;. 

A set-up, named the Low-Energy Neutrino Factory 
(LENF), has been proposed as a more conservative op- 
tion for large values of 6*13 [13, [M, [43[ . The first proposal 
for the LENF used a single baseline of 1300 km, corre- 
sponding to the Fermilab to DUSEL distance, and conse- 
quently, a lower muon energy at around 4.5 GeV [27ll38l 
|43j . A detector with good energy resolution and low en- 
ergy threshold allows the LENF to exploit the rich oscil- 
latory pattern and to achieve an important sensitivity to 
the oscillation parameters. In the original proposal, the 
detector of choice was a Totally- Active Scintillator De- 
tector (TASD) magnetized by means of a large magnetic 
cavern [38|. Other possible detectors include a magne- 
tized liquid-argon time projection chamber, which would 
be ideal due to its large size and excellent detector per- 
formance [34j . The possibility of using a non- magnetized 
detector has also been studied [iij . In this case, al- 
though the event separation between wrong- and right- 
sign muons is impossible at an individual event level, 
the use of statistical techniques allows the two chan- 
nels to be partially distinguished and good sensitivity 
can be obtained thanks to the high event numbers which 
are associated with very large non-magnetized detectors. 
Additional studies are required in order to fully under- 
stand the capability of statistical separation of right- and 
wrong-sign muons and the impact of backgrounds. Stud- 
ies of the LENF using a Magnetized Iron Neutrino Detec- 
tor (MIND) have also shown promising performance [42j 
and could serve as an intermediate step in the devel- 
opment of a future higher-energy facility. However, this 
set-up is optimized for large energies due to the poor low- 
energy efficiency of the MIND detector and consequently 



the performance typically indicates a generic preference 
for high energies. A LENF with a TASD has also been 
considered as part of an incremental neutrino factory pro- 
gram [45]. In this setting, the optimal baseline for a 
stored-muon energy of = 4.12 GeV was found to be 
around L — 1100 km. 

The conventional high-energy neutrino factory 
(HENF) and the LENF are conceptually very similar 
and are only truly separated by a choice of energy scale. 
This suggests that the two designs for the neutrino 
factory should be seen as part of a continuum in energy 
and distance. Especially in light of the recent mea- 
surement of large #13, which has led to the low-energy 
design becoming the preferred experimental option, it is 
important to perform a "green-field" study of the LENF 
in order to understand the dependence of the sensitivity 
to the oscillation parameters on the stored-muon energy 
and the choice of baseline distance and, if possible, to 
answer the question: what is the optimal LENF?. In 
this paper, we have performed this study considering 
ranges for the energy and baseline which interpolate 
between the HENF and the LENF and providing an 
understanding of how the NF design transitions between 
the two regions. In particular, we have conducted a 
finely-grained scan over the energy and baseline range 
and have identified new performance indicators for two 
primary questions: the discovery of CP-violation and of 
the mass hierarchy. The analysis that we present has, for 
the most part, been performed without fixing the value 
of #13. For the measurements of CP- violation and mass 
hierarchy this allows us to present results for the full 
range of possible values of 6*13. For completeness, we have 
also computed the discovery potential for measurements 
of 0i3 itself, as the discussion of this measurement allows 
the relevant physical effects governing the performance 
and sensitivities of the neutrino factory to be explored 
and better understood. 

After the completion of our work, a study [46| ap- 
peared which also addresses the optimal baseline and 
stored-muon energy of the LENF. Our work differs from 
Ref. in a few significant regards. First of all, we 
have performed a full, simultaneous scan over the E^-L 
parameter space and as such have a diminished reliance 
on the interpolation of results. Secondly, our study con- 
siders three different detector options which helps us to 
fully understand the potential performance of the LENF. 
Finally, in this study we have chosen 10 22 total useful 
muon decays, a range of baseline distances from 1000 km 
to 4000 km and stored-muon energies from 4 GeV to 
25 GeV. These choices were made to connect previous 
optimization studies of the LENF (1300 km, 4.5 GeV) to 
those of the HENF (short baseline - 4000 km, 25 GeV). 
Understanding this parameter space in detail allows the 
two designs to be viewed as part of a continuum and 
helps us to unite the previous work on the two designs in 
a common framework. 

This paper is organized as follows. In Section [TT] we 
discuss the details of our simulation and the assumptions 
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we have made in our model of the LENF. In Section|TTT]we 
present the results of our study and provide an analysis 
of the sensitivities of the LENF across our parameter 
space. In Section HV] we discuss the performance of the 
LENF in light of the recently measured value of #13 and 
111 Section El we summarize our findings and make a few 
concluding remarks. 



II. SIMULATION DETAILS 

In this study, we have performed simulations of the 
physics performance of the LENF over a range of experi- 
mental configurations. As our indicators of performance, 
we have computed the sensitivity of the LENF towards 
three of the most important potential discoveries for the 
future generation of neutrino oscillation experiments: the 
discovery of non-zero #13, the discovery of CP- violation 
arising from the Pontecorvo-Maki-Nakagawa-Sakata [I?} 
matrix and the determination of the neutrino mass hi- 
erarchy. We have studied the potential for discovery of 
each of these fundamental phenomena over a range of 
stored-muon energies given by 4 < < 25 GeV and 
baseline distances by 1000 < L < 4000 km. This range 
connects the regions of parameter _space traditionally as- 
sociated with the LENF desig n wl |H, HI] to those of 
the conventional HENF set-up |31j . 

It has been suggested that 1.4 x 10 21 muon decays per 
year is an attainable goal for the lower-energy accelera- 
tor facility [H[ and this estimate has been incorporated 
into previous studies of the performance of the LENF 
[43l l49j . In this work, so as to aid comparison with stud- 
ies of the standard neutrino factory [4l|, |42| . we have 
assumed 1.0 x 10 21 useful muon decays per year per po- 
larity and a run-time of 10 years which is divided evenly 
between the two polarities (10 22 useful muons in total). 
This value is in accordance with the estimates for the 
conventional neutrino factory, assuming 10 7 operational 
seconds per year [3l|. We have additionally performed 
simulations assuming the optimized value of 1.4 x 10 21 
muon decays per year per polarity, the data indicate a 
predictable uniform increase in performance and have 
been made available online [50]. In Ref. [46], the total 
number of useful muon decays was taken as 2.5 x 10 22 . 
This leads to a further increase in statistics and conse- 
quently in the discovery potential of the LENF. To un- 
derstand the performance of the LENF in the context 
of similar experimental configurations it is important to 
ensure a fair comparison is made. For this reason, in the 
analysis that follows we have assumed the conventional 
10 22 total useful muon decays. 

The detectors in our simulations measure muons aris- 
ing through both the disappearance (v^ — > or — > 
z7 M ) and the appearance channels (u e — > F M or v e — > v^). 
We have chosen not to include the platinum channel (the 
observation of v e and V e ) as previous work has shown that 
it offers only marginal improvement of the oscillation pa- 
rameter sensitivity 43]. The production of t + (t~) in 



the detector by the charged-current interactions of inci- 
dent V T {y T ) leads to the problem of tau contamination 
[HHIEij]. r-leptons have a lifetime at rest of 2.9 x 10~ 13 s 
and decay inside the detector into muons with a branch- 
ing ratio of around 17% (53[. This effect leads to an 
increased number of both wrong- and right-sign muons. 
Due to the form of the oscillation probabilities the num- 
ber of additional muons in the right-sign channel (e.g. 
from — > v T for fjT in the storage ring) are consider- 
ably larger than in the wrong-sign channel (e.g. from 
v e — > v T for /i~ in the storage ring). Subsequently, the 
effect of tau-contamination is most pronounced for mea- 
surements which rely on an accurate determination of the 
disappearance channels, for example in investigations of 
#23-m.aximality [Hl|. This contamination can lead to se- 
rious systematic uncertainties if unaccounted for. How- 
ever, it has been shown that correctly incorporating this 
additional source of muons into the analysis of the golden 
channel, via migration matrices, can resolve the system- 
atic deviations 54j. This has been confirmed in a recent 
study on the performance of the standard NF [42j where 
the change in sensitivity produced by correct incorpora- 
tion of the contamination channel was found to be small. 
The effect of tau-contamination is also expected to be 
smaller for lower-energy facilities as the number of neu- 
trinos with energies above the tau production threshold 
will be reduced. For these reasons, we have omitted the 
contamination channel in our study. 

Our simulations were performed numerically using the 
GLoBES package [55|, [56j which incorporates the Pre- 
liminary Reference Earth Model [U HH for the com- 
putation of the matter density along the baseline. In 
this package, the oscillation probabilities are computed 
via numerical diagonalization of the full Hamiltonian as- 
suming three neutrino flavors. However, it is convenient 
during the analysis of our results to introduce an ap- 
proximate expression [59j for the golden-channel oscilla- 
tion probability which is valid up to second-order in #13, 
a = Am^/Am 2 !, Am^L/E and Am^/EA 
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The first summand in this expression is referred to as 
the atmospheric term and depends quadratically on #13. 
The CP term is second and introduces dependence on 
S. The remaining part is called the solar term, which 
for small values of 6*13 can dominate the oscillation prob- 
ability. As the solar term is independent of #13, 5 and 
sign(Am§ 1 ), it can lead to a significant loss of sensitivity 
in the measurement of these parameters if dominant. 
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The optimal detector technology for the LENF has not 
yet been identified and a number of candidate designs re- 
main viable. To facilitate a comparison between the al- 
ternative designs, we have performed our simulations for 
three detectors simultaneously: a Totally-Active Scin- 
tillator Detector (TASD) and two different liquid-argon 
detectors which have optimistic and conservative perfor- 
mance estimates respectively. The TASD concept has 
been successfully implemented in the MINERVA experi- 
ment and a larger scale device has been selected for 
the upcoming NOi/ A [ 61j design. Our model of the TASD 
is based upon Ref. [38| and has a fiducial mass of 20 kton, 
35 variable- width energy bins and a constant energy reso- 
lution of 10% for both quasi-elastic and non-quasi-elastic 
events. The efficiency rises linearly from 73% to 94% over 
a range of 0.5 GeV to 1 GeV and then remains constant 
for higher energies. The background on the golden chan- 
nel is taken as a constant fraction of 1 x 10~ 3 of the events 
arising from neutral-current interactions and the same 
fraction of events from the disappearance channel which 
accounts for instances of charge misidentification. Both 
of the liquid-argon detectors are based upon parameters 
first reported in Ref. [62| and elaborated on in subse- 
quent optimization studies [43| by having a fiducial mass 
of 100 kton, an energy resolution on quasi-elastic events 
of 10% and a flat detection efficiency of 80%. The con- 
servative (optimistic) model has 22 (35) variable-width 
energy bins and an energy resolution of 20% (10%) on 
non-quasi-elastic events; the backgrounds are taken as a 
fraction of 5 x 10~ 3 (1 x 10~ 3 ) of events from both the 
neutral-current and disappearance channels. 

Although a number of large mass-scale liquid-argon de- 
tectors have been proposed 37, 63- 65] , only a few designs 
[37l [HI discuss extensions to 100 kton. The magnetiza- 
tion of large- volume detectors, of the scale considered 
in our simulations, is a particular challenge and further 
research is needed to fully assess the feasibility of the de- 
sign. Our choice of such large detector volumes (100 kton 
and 20 kton for liquid-argon and TASD respectively) is 
designed to provide an optimistic performance estimate 
which covers the full range of potentialities of the LENF 
if these technical difficulties can be overcome. It is worth 
noting that, in our approximation, the performance of the 
LENF will only depend upon the detector mass through 
its exposure (number of muon decays x fiducial detec- 
tor mass). Using this equivalence, the performance of a 
100 kton detector with 1.0 x 10 21 muon decays per year 
is expected to be comparable to a 70 kton detector with 
a larger flux of 1.4 x 10 21 decays. 

We have assumed normal hierarchy to be true through- 
out our simulations and the known oscillation param- 
eters were chosen to be sin 2 26>i2 = 0.3, 6\ 



Am 2 12 



'23 = tt/4, 

.0 x 10~ 5 eV 2 and |Am? 3 = 2.5 x 10~ 3 eV 2 . 
The uncertainty on these values was accounted for by al- 
lowing the parameters to vary during the minimization 
procedure: we allowed an uncertainty of 4% and 10% for 
the solar and the atmospheric parameters, respectively. 
These parameter choices were made in accordance with 



previous optimization studies of the LENF 43j and are 
close to the best-fit values and uncertainties from recent 
global analyses of the existing neutrino oscillation data 

la- 



in. RESULTS AND ANALYSIS 

Here we present the performance of the LENF by con- 
sidering its ability to make three major discoveries. We 
start by computing the sensitivity to non-zero #13 which, 
in spite of recent experimental results, is included to high- 
light a number of general physical effects which influence 
the performance of the the neutrino factory. After this, 
we focus on the discovery of S 6 {0, n} and Am 2 ^ > 0. 
A discussion relevant to the specific case of large 6*13 is 
presented in Section UVl 

In this work we compute the discovery potential of non- 
zero 0i3 which is defined as the exclusion of #i 3 — at 
the 3cr confidence level, marginalizing over all other os- 
cillation parameters. For the discovery of CP-violation, 
we consider the exclusion of all parameter sets with 
S E {0, 7r} (the CP-conserving values) and for the deter- 
mination of the mass hierarchy, we require the exclusion 
of all parameter sets with the incorrect hierarchy. To 
simplify the presentation of our results we use the con- 
cept of the discovery fraction, defined as the fraction of 
true values of 6 at which discovery of the related quan- 
tity is possible for a given value of #i 3 . For each choice 
of detector, baseline distance and stored-muon energy we 
have computed the discovery fraction as a function of the 
true value of 6*13 for the three quantities of interest. The 
full set of discovery plots computed in this study, show- 
ing the regions of sin 2 2613-S parameter space for which 
discovery is possible, has been made available online [50| . 

To aid our analysis, we will briefly recall some of the 
generic factors which influence the variation in perfor- 
mance of the LENF over the L-E^ parameter space. 
Maintaining a high flux of neutrinos is essential in long- 
baseline experiments as this determines the number of 
events that can be observed at the detector and governs 
the statistical significance of any observations. The total 
flux incident on the detector increases with the stored- 
muon energy as E 2 , whilst decreasing with longer base- 
lines as L~ 2 . The observed number of wrong-sign muons 
is heavily influenced by the oscillation probability, Eq. [TJ 
which introduces an additional dependence on the base- 
line distance and on the energy of the individual neu- 
trinos, E v . Observing events which come from the first 
oscillation maximum is important as it ensures a large 
signal in the appearance channel. Events from the low- 
energy part of the spectrum contain important informa- 
tion on CP-violation, as for these values of E v the os- 
cillation probability exhibits a strong dependence on the 
CP-violating phase 5. Matter effects lead to an enhance- 
ment or suppression of the oscillation signal compared 
to the same process in vacuum. Observing the sign of 
this difference can provide us with information on the 
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FIG. 1. The #13 discovery range as a function of baseline, L, and stored-muon energy, E^. The left-hand column shows the 
lowest value of sin 2 20n for which the discovery fraction is 100% for all higher values and the right-hand column shows the 
lowest value of sin 2 2#i3 for which there is a non-zero discovery fraction. The intensity of each point is given by the value of 
log 10 (sin 2 20i3 ). 
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neutrino mass hierarchy. In Eq. [T] matter effects are 
present if A ^ and their influence increases with base- 
line distance and neutrino energy. Sensitivity to a given 
measurement can also be affected by the presence of de- 
generacies, the locations of which generally depend on 
both L and E v . At the detector, the neutrino-nucleon 
interaction cross-sections depend on the energies of the 
incident neutrinos. Flavor-tagged detection is made pos- 
sible for charged-current interactions by observing the 
emitted charged lepton. The total cross-section for this 
channel increases with energy [66l.l67| and, assuming that 
all else is kept equal, leads to an improved number of 
events at higher-energy facilities. 

To understand the variation in performance of the 
LENF as we vary L and E^, the individual dependences 
mentioned above must be considered in combination. In 
the following sections, we will discuss how these effects 
can explain the sensitivity of the LENF to non-zero #13, 
CP violation and the mass hierarchy. 



A. Discovery of non-zero #13 

Here we present the results of our simulation which ad- 
dress the ability of the LENF to discover non-zero #13. In 
light of the recent measurements of #13 , optimization to- 
wards this measurement has become largely peripheral. 
However, the results in this section, and the accompa- 
nying discussion, introduce a number of elements which 
govern the performance of the LENF in general. As such, 
it is instructive to first understand this measurement be- 
fore moving on to the question of CP-violation discovery 
and the determination of the mass hierarchy. A selection 
of our results can be seen in Fig.[T]where the left column 
shows, as a function of L and E^, the lowest value of 
sin 2 2#i3 for which discovery of #13 7^ can be expected, 
independently of S, for all higher values. This condi- 
tion is equivalent to ensuring that all higher values of 
sin 2 20i3 have discovery fractions of 100%. On the right, 
Fig. Q] shows the lowest value of sin 2 2#i3 which can be 
discovered for at least one value of 5, which corresponds 
to a 0% discovery fraction. These can also be seen as 
the end points of the #13 discovery fraction curves seen 
in previous studies (see for example Fig. 2 in Ref. [42|). 
Physically, the first of these quantities gives the small- 
est value of #13 above which we expect a discovery and 
provides a conservative estimator of performance. The 
second quantity is a complementary optimistic estima- 
tor and tells us the smallest value of #13 at which we 
could possibly make a discovery. The range of these two 
parameters gives the region of intermediate performance 
where we find discovery fractions between zero and one; 
discovery in this region is dependent on the exact value 
of 5. 

In Fig. Q] we see that the TASD expects #i3-discovery 
to at least sin 2 29^ > 10~ 3 and has the possibility of 
extending this limit by an order of magnitude. In com- 
parison, the optimistic liquid-argon detector can discover 



non-zero $13 down to at least sin 2 2#i3 > 10 4 and pos- 
sibly as low as sin 2 26»i 3 > 3 x 10~ 5 . The TASD and the 
conservative liquid-argon detector generally offer compa- 
rable sensitivities, which are both worse than those of 
the optimistic liquid-argon detector. For each detector, 
we attribute this to different causes. The conservative 
liquid-argon detector sees a similar total number of events 
as the optimistic liquid-argon detector however its higher 
backgrounds lead to a poor signal to background ratio 
which reduces its sensitivity. In contrast, the TASD has a 
superior control of the dominant backgrounds to the con- 
servative liquid-argon detector but sees far fewer events 
due to its smaller size. These two effects reduce the at- 
tainable experimental sensitivity by a similar degree. 

To further the analysis of our results, we will discuss 
the four quarters of this parameter space separately. In 
the short-baseline and low-energy region (SB-LE) with 
L < 2500 km and E^ < 14 GeV, we see 100% discov- 
ery fractions for 10~ 4 < sin 2 2#i3 < 10~ 25 depending 
on the choice of detector technology. In this region, the 
atmospheric term is relatively suppressed in the P eM os- 
cillation probability. This suppression is a result of small 
6*i3 and the relative enhancement of the CP and solar 
terms at lower neutrino energies. This leads to a poor 
performance as the signal becomes decreasingly sensitive 
to #13 and the CP term introduces a more complicated 
dependence on the oscillation parameters. The decline in 
performance found towards the very shortest baselines is 
a consequence of the reduction in oscillation probability 
at small L. This leads to a poor signal-to-background 
ratio as the number of wrong-sign muons decreases. In 
the region of short baselines and high energies (SB-HE) 
with L < 2500 km and E^ > 14 GeV, the larger ener- 
gies increase the relative importance of the atmospheric 
term and we see that the situation is marginally improved 
with respect to the SB-LE region with 100% discovery 
fractions for 10~ 4 < sin 2 2#i3 < 10 -3 depending on de- 
tector technology. The overall improvement in this region 
can be seen as simply a leveling-off of the higher energy 
improvements of the SB-LE regions: for a fixed base- 
line distance, there is negligible improvement to be found 
when moving from E^ — 15 GeV to E^ — 25 GeV. This 
plateauing effect is associated with the energy spectrum 
of the neutrino factory, which rises almost linearly from 
E u = to its peak at E v k, 2E fl /3 and then drops sharply 
at the kinematic cut-off E v = E^. As we increase the 
stored-muon energy, the approximately linear tail of this 
spectrum only decreases slightly and therefore increasing 
the stored-muon energy at a LENF can be thought of 
as simply adding higher-energy neutrino events on top 
of the previous low-energy spectrum. Consequently, the 
difference in performance between an experiment with a 
low and a high stored-muon energy can be estimated by 
determining the importance of the neutrino events oc- 
curring with energies between the two stored-muon ener- 
gies. Moving up through the SB-HE region, the informa- 
tion provided by the additional high-energy neutrinos is 
decreasingly useful because for higher-energy neutrinos, 
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the L/E ratio is smaller and the oscillation probability 
decreases. This effect leads to a law of diminishing re- 
turns, where the discovery reach remains approximately 
constant. However, as we move to higher energies, the 
signal to background ratio decreases and with the nec- 
essary introduction of additional backgrounds and the 
slight decrease in the low-energy part of the spectrum, 
the sensitivity in this region is expected to ultimately be 
reduced. 

For the regions of parameter space with baselines of 
L > 2500 km, we again divide the parameter space 
along = 14 GeV into two quarters: the long-baseline, 
low-energy region (LB-LE) and the long-baseline, high- 
energy region (LB-HE). The LB-HE region has in gen- 
eral the best 100% discovery fractions of all of the pa- 
rameter space. In this region, the design is approaching 
the HENF configuration where the signal to background 
ratio is increased and, thanks to the energy-dependent 
relative suppression of the solar and CP terms, the at- 
mospheric contribution to the oscillation probability has 
a significant influence on the appearance channel. This 
permits the true value of #13 to be smaller than in the 
SB-LE region whilst still providing an appreciable sig- 
nal and therefore furthering the discovery reach. In con- 
trast, the LB-LE region displays the poorest performance 
across the parameter space. We find 100% discovery frac- 
tions at around sin 2 20 i3 w 10~ 25 for the TASD and 
conservative liquid-argon detectors, whilst the optimistic 
liquid-argon detector has a 100% discovery fraction for 
10 -3 < sin 2 2#i3 % 10~ 4 . The poor performance in this 
region can be explained as an effect of low statistics: the 
neutrino flux is reduced as the baseline increases due to 
the effects of dispersion on the beam. In the LB-LE re- 
gion this effect is compounded with a small neutrino flux 
at production. Consequently, the number of wrong-sign 
muons incident on the detector becomes increasingly re- 
stricted for these parameter choices. 

In general, the optimal configuration depends upon the 
magnitude of the true value of 6*i 3 . For large values of 
6>i3, as recently confirmed by Daya Bay and RENO[9l[l0|. 
all configurations are equally capable of confirming this 
effect. 

The 6*i3 discovery reaches reported here are consistent 
with studies performed on the conventional NF [42[ and, 
although making a direct quantitative comparison is dif- 
ficult, the discovery reaches are of a similar magnitude. 
This behavior is also reported by Ref. [46| although, due 
to a difference in exposure, once again the results differ 
quantitatively. 



B. Discovery of CP-violation 

In Fig. [5] we present the CP-violation discovery frac- 
tion as a function of baseline and stored-muon energy at 
different true values of sin 2 26>i3 for the TASD and the 
optimistic liquid-argon detector. The performance of the 
conservative liquid-argon detector is similar to that of the 



TASD and has been omitted from this paper (however, 
these plots may still be found online [50]). In the SB-LE 
region, for sin 2 26>i3 > 10~ 2 we generally see very strong 
performance with discovery fractions of 70% to 90% de- 
pending on the choice of detector technology. These sen- 
sitivities confirm the expectations behind the original 
motivation for the LENF: a facility with the ability to 
measure low-energy neutrino events has access to the os- 
cillation spectrum near the second maximum where CP- 
violating effects are most pronounced. The simulation 
has also been performed for smaller values of sin 2 26*13 
where the performance of the SB-LE region starts to 
worsen as the (^-dependent terms in the probability are 
increasingly suppressed. At the lowest energies, this sup- 
pression is compounded by the enhancement of the #13- 
and ^-independent solar term. For sin 2 26*13 w 10~ 4 there 
is only negligible coverage for the TASD however the op- 
timistic liquid-argon detector still maintains a discovery 
fraction of around 60% due to its combination of high 
statistics and strong signal to background ratio. Consid- 
ering experiments with higher stored-muon energies, we 
see scant improvement in the SB-HE region as the ad- 
ditional events at high energy provide little information 
on the parts of the oscillation spectrum which exhibit 
the most sensitivity to the CP-violating phase. Gener- 
ally we see discovery fractions of around 70% to 90% 
for sin 2 26»i 3 > 10~ 2 , which drops to 60% to 70% for 
sin 2 26*13 ~ 10~ 3 . For the smallest values of #13, the 
TASD has once again negligible sensitivity whilst the op- 
timistic liquid- argon detector can still determine the ef- 
fect of CP-violation in 60% to 70% of cases. 

As in Section IIII Al while the baseline distance is in- 
creased there is a decrease in event numbers due to a 
weakening of the neutrino flux arising from long baselines 
and this may be compounded by an additional weakening 
of the flux at low energies. We generally find the lowest 
discovery fractions in the LB-LE region and this contrast 
is especially marked in the case of small #13. However, 
for long baselines but high energies (LB-HE) we see good 
sensitivity to CP-violation, especially when #13 is small. 
In this region, both the neutrino flux at production and 
the neutrino-nucleon cross-sections are increased and this 
helps to mitigate the effect of baseline distance on the 
event numbers. The additional influence of appreciable 
matter effects over longer baselines and the inclusion of 
neutrinos which probe the most CP-sensitive parts of the 
oscillation spectrum further improve the sensitivity. For 
sin 2 26*13 ^ 10~ 3 , the LB-HE region has comparable dis- 
covery fractions to those of the SB-LE region however 
for sin 2 26*13 w 10~ 4 the only significant sensitivity is to 
be found in the LB-HE region with discovery fractions 
of 40% to 70% for the TASD and optimistic liquid-argon 
detectors, respectively. 

As mentioned previously, the true dependence of our 
simulations on detector mass, runtime and the number 
of useful muon decays per year is through their product, 
referred to as the exposure. In Fig. [3J we can further our 
understanding of the CP-violation discovery fraction by 
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FIG. 2. CP- violation discovery fractions as a function of baseline, L, and stored-muon energy, _E M . Each column (row) shows 
the discovery fraction for a different 'true' value of sin 2 2#i3 (detector option) which is given at the top-left (top-right) of each 
plot. 



considering its dependence on this parameter. In gen- 
eral, for large values of #13 the discovery fraction reaches 
a plateau for each experimental configuration. This limit- 
ing value represents the inherent limitations of the exper- 
iment: exposure directly determines the statistics of our 
experiment and at some point the measurement uncer- 
tainties will be dominated by systematic effects for which 
an increase in statistics can confer only slight improve- 
ment to the experiment's sensitivity. The plots in Fig. [3] 
show the difference in performance for a modest (10%, 
left plot) and a more optimistic (1%, right plot) energy 
resolution over a range of values of a uniform systematic 
error on the signal and backgrounds. We see that vari- 
ations in magnitude of the systematic errors induce the 



greatest change in the attainable CP-violation discovery 
fraction. These effects can lead to a significant decline 
in performance; for 10% (1%) energy resolution, there is 
a decrease in the discovery fraction of 8% (6%). This 
limiting influence on the discovery fraction which arises 
through systematic uncertainties is quite stable to vari- 
ations under baseline and energy, assuming that these 
choices do not generate a significant probabilistic sup- 
pression of the number of events. This effect helps to 
explain the observed uniformity in parts of Fig. [21 for 
example in the SB-HE region. We have also investigated 
the impact of alternative sources of systematic limita- 
tions. For instance, reducing the prior uncertainty associ- 
ated with the matter density can lead to modest improve- 
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FIG. 3. The CP-violation discovery fraction as a function of exposure for a range of systematic errors. The plot on the left 
(right) shows the discovery fraction assuming an energy resolution of 10% (1%). The scale of the horizontal axis is chosen to 
coincide with runtime assuming a 20 kton detector and 10 21 useful muon decays per year per polarity. These plots use the 
TASD and have #13 set to the Daya Bay best-fit value of sin 2 2^13 = 0.092 All otherwise unspecified parameters are as 
described in Section [TT1 



merits in sensitivity; however, we find that the unilateral 
improvement of any one systematic factor leads to little 
impact on the sensitivity obtained in Fig. [3j it is neces- 
sary to reduce all systematics uniformly to significantly 
further the physics reach. Generally for the parameter 
ranges that we have studied, it is the energy resolution 
and overall systematic error that are responsible for the 
greatest variation in the attainable CP-discovery frac- 
tion. 

In common with the discovery of #13 ^ 0, the opti- 
mal configuration of the LENF for CP-violation discov- 
ery divides into two scenarios depending on the size of 
6*i 3 . For sin 2 28\ 3 > 10~ 3 , provided the LB-LE region is 
avoided, the CP-violation discovery fractions are almost 
uniform and the exact configuration is unimportant. As 
#13 decreases, the (^-dependent signal is suppressed and 
appreciable sensitivities can only be found in the LB-HE 
region where the advantages of the HENF design start to 
become relevant. 

The per formance of the single-baseline HENF has been 
shown [42] to share generic features with our data: the 
SB-LE region suffers from lower discovery fractions and 
once #13 has decreased to around sin 2 26*13 = 10 -4 , the 
only sensitivity can be found in the LB-HE region. The 
CP-violation discovery fractions of the single-baseline 
HENF are very similar to those of the LENF with the 
TASD whereas, for all values of sin 2 2613 that we have 
studied, the liquid-argon detector has discovery frac- 
tions higher by around 10%. For a two-baseline HENF, 
with the longer baseline chosen to be at the magic base- 
line [42j, the comparison changes depending on the size 
of sin 2 26>i 3 . For sin 2 20 13 > 10" 2 the HENF performs 



similarly to the LENF with TASD whilst the optimistic 
liquid-argon detector has higher discovery fractions, once 
again by around 10%. For values of #13 in the range 
10" 4 < sin 2 26»i3 < 10" 3 , the HENF starts to out- 
perform the LENF with TASD. However, the perfor- 
mance of the LENF with an optimistic liquid-argon de- 
tector remains comparable. Our results in the SB-LE re- 
gion agree qualitatively with those computed in Ref. [43[ 
for a LENF with L = 1300 km and = 4.5 GeV and 
also agree qualitatively with the recent simulations of the 
LENF HU. 



C. Determination of the Mass Hierarchy 

In order to present the results of our simulations re- 
garding the determination of the mass hierarchy, we have 
computed analogous quantities to those used in Sec- 
tion IIII A[ In Fig. 2] we have plotted, as a function of 
L and E^, the smallest value of sin 2 26*13 at which hier- 
archy determination is expected to be possible: the 0% 
discovery reach. In Fig. @] we show the smallest value 
of sin 2 2#i3 for which the mass hierarchy can be deter- 
mined for all higher values of sin 2 2#i3 , independently 
of the value of 5: the 100% discovery reach. We see 
that, for most of the parameter space, the TASD pro- 
vides 100% discovery reaches of around sin 2 26*i 3 > 10~ 2 
and 0% discovery reaches that are smaller by an order 
of magnitude. The optimistic liquid-argon detector of- 
fers a similar 100% discovery reach for baselines below 
L = 2500 km but, for baselines greater than this, can pro- 
duce a significantly lower limit between sin 2 26*13 ^ 10~ 3 
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FIG. 4. Hierarchy determination as a function of baseline, L, and stored muon energy, for the TASD (top row) and 
the liquid- Argon detector with optimistic performance estimates (bottom row). The left column shows the most conservative 
discovery limit: the lowest value of sin 2 2#i3 for which all higher values allow determination of the mass hierarchy independently 
of 8. The right column shows the optimistic limit: the lowest value of sin 2 2#i3 for which there is a non-zero discovery fraction 
for the exclusion of the incorrect hierarchy. In each plot, the black lines describe regions for which discovery is possible up to 
the value of sin 2 2#i3 as labeled. The intensity at each point is given by log 10 (sin 2 2^13) of the relevant discovery limit. 



and sin 2 26>i3 > 10~ 4 . The 0% discovery reach in this 
case is generally sin 2 26>i3 > 3.2 x 10~ 4 going down to 
sin 2 26*13 ^ 3.2 x 10~ 5 for the longest baselines at around 
L = 3500 km. It is clear that these plots exhibit a 
stronger dependence on baseline distance than has been 
seen in previous plots. This was to be expected as it 
is well known that matter effects are crucial in lifting 
the hierarchy degeneracy and that these effects increase 
with longer baselines. This can be seen in the golden- 
channel oscillation probability by considering the differ- 
ence between the appearance probability for wrong- and 



right-sign muons. Depending on the neutrino mass hier- 
archy, one of these channels is suppressed and the other 
enhanced: this discrepancy grows with longer baselines. 
An exception to this pattern is found in the LB-LE re- 
gion where there is a notably poor performance compared 
to the other regions: the particularly low neutrino flux 
arising from the combination of long baselines and low 
stored-muon energies leads to this decreased sensitivity. 
These results are in qualitative agreement with predic- 
tions for the standard NF [42| and with previous studies 
of the LENF [Hill. 
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It has been suggested _68] that a low-energy neutrino 
factory with a "bimagic" baseline of around L = 2540 km 
and stored-muon energy of = 5 GeV would offer 
a pronounced sensitivity to the neutrino mass hierar- 
chy. This claim was motivated by studies of superbeams 
[63 . [rol l which looked for points of L-E^ parameter space 
which show a clean separation between the hierarchy- 
conjugate oscillation probabilities. Such a point was 
found at L — 2540 km and E = 3.3 GeV where the oscil- 
lation probability for the inverted hierarchy is small and 
independent of 5. In contrast, the oscillation probability 
for normal hierarchy is much larger and consequently the 
neutrino fluxes for different mass hierarchies arc expected 
to differ significantly at this point. The same property is 
also found at L — 2540 km and E = 1.9 GeV but this 
time the probability for the normal hierarchy is small and 
(^-invariant whilst that of the inverted hierarchy is large. 
It was the existence of two "magic" energies for the same 
value of L that lead to the idea of a "bimagic" baseline. 
A neutrino factory with a stored-muon energy of 5 GeV 
and a baseline of L = 2540 km would produce a spectrum 
covering both of these points and this has been shown to 
lead to a strong sensitivity to the neutrino hierarchy [68j . 
In Fig. [5] we show our results indicating how the discov- 
ery reach depends on baseline distance for a selection of 
muon energies similar to the bimagic set-up. At a muon 
storage energy of E^ = 5 GeV we see evidence for a 
minimum in the 0% discovery limit at baseline distances 
around 2500 to 2800 km. This lower limit corresponds 
to the optimistic performance of the LENF and it is im- 
portant to note that there is no corresponding minimum 
in the conservative estimate shown by the the 100% dis- 
covery limit. We see from the other bands in Fig. [5] how 
this feature changes as the stored-muon energy increases: 
the minimum flattens out and drifts to higher baselines. 
We see that, for all stored-muon energies, a facility with a 
baseline distance below L = 2500 km can improve its dis- 
covery reach notably by increasing its baseline to at least 
2500 km. Beyond this, it appears that if muon storage 
energies higher than 5 GeV are available then the bimagic 
choice is not the optimal configuration as the discovery 
reach can be additionally furthered by increasing both 
the energy and baseline. Furthermore, it is important to 
remember that even if stored muon energies are fixed at 
5 GeV, a baseline of L w 2500 km only maximizes the 
optimistic performance of the design and any potential 
sensitivity advantage would crucially depend on the value 
of S. 

The behavior in Fig. [5] can be understood in light of 
our previous analysis. Matter effects are necessary to 
lift the hierarchy degeneracy and these are increased sig- 
nificantly by the use of long baselines. The reason that 
this trend does not appear for the E^ = 5 GeV case is 
because, at such a low energy, the experiment becomes 
statistically limited in the LB-LE region and necessitates 
the move to higher energies. The fact that the configu- 
rations near the bimagic baseline do not confer equal im- 
provements for the conservative, 100% discovery limit is 



a result of the construction of the bimagic baseline crite- 
ria: the degree of contrast between hierarchy-conjugate 
probabilities is not ensured to be large and its magni- 
tude depends on the exact value of S. Although for some 
values the contrast is sizable, for others it is greatly re- 
duced. For example, at the bimagic configuration, when 
<5 = the lowest magic energy at E v — 1.9 GeV predicts 
rates which are almost identical for both hierarchies and 
can offer little discriminatory information. The presence 
of values of 5 for which the hierarchy distinction is less 
marked leads to a weaker 100% discovery limit, which is 
designed to measure exactly this worst case scenario. 

Provided the LB-LE region is avoided, identifying the 
optimal configuration for measurements of the mass hi- 
erarchy reduces to the observation of a simple correla- 
tion between the expected magnitude of sin 2 2#i3 and 
the baseline distance. For the smallest values of #13, con- 
figurations in the LB-HE region are necessary. However, 
for the larger values recently measured by Daya Bay and 
RENO, the physics reach is quite stable and the exact 
configuration in our parameter space is unimportant. 



IV. DISCUSSION 

In light of the latest experimental results [6-10], the 
optimization of the LENF for large values of #13 has be- 
come essential. The Daya Bay and RENO experiments 
can both independently exclude #13 — at around 5ct and 
suggest an allowed region of around 0.07 < sin 2 2#i3 < 
0.14 at la [9j, llOj. Our simulations show that these val- 
ues lie in the optimal region for the performance of the 
LENF: CP-violation discovery fractions of 80-90% are 
attainable in most of the parameter space with an opti- 
mistic liquid-argon detector or TASD. This sensitivity 
can be understood by the low-energy enhancement of 
the CP term in the golden-channel oscillation probabil- 
ity combined with a relative suppression of the solar term 
because of the large value of #13. Additionally, the mass 
hierarchy will be measurable independently of the choice 
of 5 for all of the detectors, baselines and stored-muon 
energies that we have considered in this study. This capa- 
bility arises because of the effect of the matter potential 
which generates an enhancement or suppression in the 
expected number of wrong-sign muons depending on the 
neutrino hierarchy. For a large value of $13 , a significant 
number of oscillated events should occur and it is un- 
likely that the parameter sensitivity will be statistically 
limited. As the motivation for detectors with very large 
masses is primarily one of statistics, a large value of 613 
allows a reduction in detector mass whilst maintaining a 
similar level of performance. Although statistical errors 
are less relevant with a larger mixing angle, systematic 
errors become increasingly important. A careful study 
of the systematic uncertainties will be an important next 
step in the oscillation analysis. 

For a related discussion of the effect of large values of 
013 on the sensitivities of the LENF, see also Ref. [461 ]. 
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V. CONCLUSIONS 

The Neutrino Factory has previously been studied in 
two distinct configurations: the conventional HENF uses 
25 GeV muons and has two baselines at 4000 km and 
7000 km. More recently, the idea of the LENF has been 
presented: this uses low-energy muons, typically around 
5-10 GeV, and a unique baseline at 1500-2000 km. This 
set-up exploits the rich oscillatory pattern of the appear- 
ance probability and, by mitigating the effect of degen- 
eracies among the unknown neutrino parameters, can 
provide an excellent physics reach for larger values of 
(9 13 . Motivated by recent experimental results, we have 
performed a "green-field" study of the dependence of 
the performance of the LENF on the choice of stored- 
muon energy, 4 < < 25 GeV, and baseline distance, 
1000 < L < 4000 km, in order to ultimately identify the 
configuration of the optimal LENF. 

In this paper we have presented the results of numerical 
simulations on the ability of the LENF to answer three 
questions: is #13 non-zero, does the neutrino mixing ma- 
trix give rise to CP-violating phenomena and what is the 
correct neutrino mass hierarchy? We find that non-zero 
6*i3 can be discovered for almost all values of baseline 
and stored-muon energy down to sin 2 26*13 > 10~ 3 and 



we expect limited sensitivity for a further order of mag- 
nitude. For CP-violation, discovery fractions of 70% to 
90% are expected to be attainable for sin 2 26>i3 > 10~ 2 
and of 50% to 70% for sin 2 28 13 > 10~ 3 . We have 
shown that, provided that the extremal configurations of 
our E^-L parameter space are avoided, specifically very 
short baselines or long baselines paired with low energies, 
these performance estimates are quite general and do 
not require significant fine-tuning of the baseline or the 
stored-muon energy. The neutrino mass hierarchy is ex- 
pected to be accessible to the LENF for sin 2 2#i3 > 10 -2 
with the possibility of discovery generally extending to 
sin 2 26*13 ^ 4 x 10~ 4 . For the hierarchy determination, 
the dependence on baseline is clearly seen and if base- 
lines of L > 3000 km are selected, discovery reach could 
be extended even further to sin 2 26*13 >4x 10 -5 . 



We have also considered the potential advantages of 
the bimagic baseline configuration given by L = 2540 km 
and E^ — 5 GeV. We have shown that for this choice of 
parameters, although a performance maximum is present 
in the 0% discovery reach, there is no corresponding max- 
imum in the 100% discovery reach. This means that, if 
considered conservatively, such a configuration is of lim- 
ited benefit for a NF. We have also shown that higher 
energy configurations can provide an improved discov- 
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ery reach and we see a drift in the optimal configuration 
towards higher energies and baselines. 

Now that #13 has been measured to be relatively large 
and sin 2 2#i3 ps 0.09, we see that the sensitivity of the 
LENF is comparable for the majority of baseline distance 
and stored-muon energy arrangements. Furthermore, 
we expect the results not to be statistically limited, 
implying that smaller detectors might be considered to 
provide the required sensitivities and that a detailed 
study of the impact of systematic errors is required. 
For these values of #13, the optimal LENF for the 
resolution of the most compelling open questions of 
long-baseline physics has a very broad design. This 
stability of performance shows that the LENF concept 
is versatile, with flexibility to accommodate additional 
design criteria, whilst also providing a strong sensitivity 
to the most important physical quantities. 
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